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Protonation of DNA basepairs is a reversible phenomenon which can be controlled by tuning the 
pH of the system. Under mild acidic conditions, the hydrogen bonding pattern of the DNA basepairs 
undergoes a change. We study the effect of protonation on the electronic properties of the DNA 
basepairs to probe for possible molecular electronics applications. We find that, under mild acidic 
pH conditions, the A:T basepair shows excellent rectification behaviour which is, however, absent 
in the G:C basepair. The mechanism of rectification has been discussed using a simple chemical 
potential model. We also consider the non-canonical A:A basepair and find that it can be used as 
efficient pH dependent molecular switch. The switching action in A:A basepair is explained in the 
light of 7T-7T interactions which lead to efficient derealization over the entire basepair. 



The research on electron transport through molecules 
is gaining a lot of attention in recent years due to 
their potential applications in a large class of nanoscale 
technologies pQ . With advances in nanotechnology, there 
has been an active effort to use biomolecules for molecu- 
lar electronic applications. This has been especially true 
for DNA based electronics owing to the unique geometry 
and stability of DNA which is built up on weak inter- 
molecular 7r-stacking and specific hydrogen bonding (H- 
bonding) interactions between Adenine: Thymine (A:T) 
and Guanine :Cytosine (G:C) basepairs[2 . 

So far, the studies on DNA conductivity have been 
mostly restricted to understanding the effect of tt-tt in- 
teractions between the DNA bases. This is due to the 
fact that experimental work on DNA conductance is car- 
ried out with DNA strands and the measurements are 
always between end to end of the DNA strand, or, in 
other words, about the helical axis of the DNA strand |3J. 
The importance of H-Bonding in these experiments have 
been restricted to the formation of the DNA double helix, 
which stems from the specific basepairing pattern among 
the DNA bases. 

In this work, we have analyzed the role played by the 
H-bonding between the DNA bases under the influence of 
mild acidic conditions. It is known that the protonation 
of the DNA basepairs under mild acidic conditions is a re- 
versible phenomenon and that it disrupts the H-Bonding 
profile between the DNA bases [4]. We find that, this 
protonation under mild conditions changes the electronic 
properties of the DNA basepairs significantly and the ef- 
fect being different for different basepairs. The effects 
of pH dependent reversible protonation-deprotonation on 
the electronic properties of a molecular system have been 
recently studied experimentally for di-block molecular 
diodes, where the inversion of the rectifying effect was 
found on selective protonation and deprotonation [5]. 

We have used the SIESTA and Gaussian 03 codes 
for our first principles DFT simulations [6j |7]. Within 
SIESTA, we have used the generalized gradient approx- 



imation for the exchange-correlation energy functional 
in the version of Perdew, Burke, and Ernzerhof [8] . A 
double-^ basis set with the polarization orbitals has been 
included for all the atoms [5]. To model the gold elec- 
trodes, we placed two 27 atoms gold clusters at either 
end of the DNA basepairs. The coordinates of the gold 
clusters were obtained by relaxing two Au(lll) surfaces 
separated by 30 A, in a periodically repeated cell. Each 
surface was comprised of three Au(lll) planes in a 2x2 
surface unit cell along the [111] direction. The clusters 
thus derived from the relaxed surface calculations were 
then attached to the DNA basepairs via a thiol linker 
(-CH 2 S). The Sulphur head group of the thiol linker was 
aligned on the fee hollow site of the Au(lll) surface us- 
ing S-Au distance of 2.51 A, previously reported in lit- 
erature [TU]. For the subsequent SIESTA calculations, 
the initial geometries of the basepairs were generated by 
Gaussian 03 calculations. Within Gaussian 03, we have 
used the polarizable continuum model with water dielec- 
tric to model the solvation effects. The large dielectric 
continuum of water reduces the electrostatic repulsions 
and stabilizes the formal positive charges introduced due 
to protonation. We have used the LANL2DZ basis set 
for obtaining the structures of the basepairs constrained 
between the two gold electrodes, wherein the gold elec- 
trodes arc modeled as an equilateral triangle with an 
Au-Au distance of 2.9 A which was obtained from the 
relaxed gold surfaces and which matches well with the 
previously reported experimental results [11 . All the cal- 
culations in SIESTA for the DNA basepairs attached to 
the gold clusters have been performed in a 30 x 30 x 60 A 3 
supercell, which is large enough for interactions between 
the neighboring fragments to be negligible. To model the 
reversible denaturation of the DNA basepairs, which can 
be controlled by pH variation, we have considered only 
the protonation of the amino nitrogen's, which can be 
easily achieved under mild acidic conditions |12j. 

Using the above described methodology, we have 
studied three different basepairs, namely, A:T (Ade- 
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FIG. 1: (color online) (a) Probable scheme for experimental 
realization. Optimized geometries of the DNA basepairs: (b) 
A:Tjvp, A:T P , (c) G:C np, G:C p , and (d) A:Ajv P , A:A P . 
The Gold clusters have been omitted for clarity. 



nine: Thymine), G:C ( Guanine :Cytosine), and A: A (Ade- 
nine: Adenine). The first two are the naturally occurring 
watson-crick basepairs where the number of H-bonds be- 
tween the bases are two for A:T and three for G:C. As a 
first observation we note that on protonation both Ade- 
nine :Thymine (A:T) and Guanine :Cytosine (G:C) base- 
pairs lose the N-H...N H-bond, however, they still retain 
the N-H...0 H-bond. Thus, for a complete loss of H- 
bonding under mild acidic conditions, we consider the 
non-canonical Adenine: Adenine (A:A) basepair, which 
loses both the N-H...N H-bonds on protonation. Here, 
we would like to point out that the systems studied by 
us can be realized by experimental synthesis. We give 
a probable scheme for the experimental realization in 
Fig. 1(a) based on the pioneering work of Mirkin et al. 
and Alivisatos et al. [T2]- The first step in the exper- 
imental method involves attaching DNA bases capped 
with thiol groups, which bind to gold, to a gold surface. 
Then the complementary bases arc similarly attached to 
the surface of a gold nanocluster via thiol linkages. The 
bases attached to the gold nanoclustres are then immo- 
bilized on the gold surface by spotting solutions. The 
immobilization is governed by the complementary base 
pairing. For the synthesis of the non-canonical basepair 
A:A, one needs to functionalize both the gold surface and 
gold cluster with Adenine. Hereafter, we denote the non- 
protonated basepairs as A:Tj\rp, G:C np, A:Anp and the 
protonated basepairs as A:Tp, G:C p, A:Ap. 

Upon optimization of the nonprotonated and the pro- 
tonated basepairs, we find that the basepairs do not get 



aligned perpendicular to the surface of the gold cluster. 
The tilt with respect to the surface is about 106° and this 
is due to the presence of the thiol linker which are used 
to attach the basepairs to the gold cluster surface. In 
Fig. 1(b), (c) and (d), we present the optimized geome- 
tries of A:TjVp,p, G:Cnp,p and A:Anp,p, respectively. 
We find that the A:Tatp, G:Catp and A:Ajvp form pla- 
nar structures due to almost linear H-bonding. However, 
the planarity of the H-bonds are lost on protonation in 
A:T P , G:C P and A:A P . 

For A:Tp (Fig. 1(b)), we find that the protonation leads 
to the formation of a bidentate H-bond between two N-H 
and one O (1.82 A, 145° and 1.92 A, 141°) as compared 
to the linear N-H...0 (1.72 A, 173°) and N...H-N (1.67 A, 
178°) H-bonds in A:T W p. For G:C P (Fig. 1(c)), the pro- 
tonation disrupts the central N-H...N H-bond and leads 
to the reorientation of the bases to form two N-H...0 
(2.02 A, 174° and 1.67A, 157°) H-bonds as compared to 
the two N-H...0 (1.86 A, 179° and 1.62A, 179°) H-bonds 
and one N-H..N (1.84 A, 178°) H-bond in G:C NP . The 
case of A:Ap (Fig. 1(d)) is distinct. The protonation in 
this case disrupts the two N-H...N H-bonds and the re- 
sultant structure is stabilized by partial tt-tt interactions 
between the Adenine bases. In A:Aatp, there are two 
equivalent N-H...N (1.84 A, 178°) H-bonds. Thus, the 
main point of our findings is that the protonation affects 
the nature and the strength of the H-bonds which leads 
to a considerable change in the structure of the basepairs 
with little change in energy as the N-H...N H-bond has 
the lowest H-Bond energy. We investigate the transport 
phenomenon of the non-protonated and the protonated 
basepairs to analyze the role played by H-bonds and the 
it-it interactions affecting the electronic structure of the 
basepairs thereby the related conductivity characteris- 
tics. 

To understand the transport properties of the systems 
under study, we have calculated the transmission func- 
tion T(E) at different bias voltages using the Green's 
function methodology fT3]. The Hamiltonian (H) and the 
overlap (S) matrices arc obtained from the self-consistent 
DFT calculations. To model the effect of bulk electrodes, 
we modify the Hamiltonian of the system by adding the 
imaginary self energy terms, Xl,p, (L and R being left 
and right electrode indices respectively) arising due to 
the bulk gold, to the orbitals of the gold cluster. Thus, 
the original Hamiltonian H gets modified to H as 



H = H 



(1) 



With the modified Hamiltonian, we calculate the Green's 
function G(E) and the transmission T(E) is then calcu- 
lated as 



T(E) = Tr[T L G{E)T R G+{E)] 



(2) 



where T^.p are the anti hermetian parts of the self-energy 
matrices which describe the broadening of the energy lev- 
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FIG. 2: (Color online) Transmission T(E) as function of ap- 
plied bias voltage, (a) G:Cjvp, (b) G:Cp, (c) A:Tjvp, (d) 
A:T P , (e) A:Ajvp and (f) A:A P . 

els of the system due to the coupling to the electrodes. 

r ijB (£0=*(E i> j l -Ei >fl ) (3) 

In Fig. 2, we plot the transmission T(E) as a function 
of bias voltage for G:C(jvp,p), A:T(jvp,p) and A:A(jvp,p)- 
We restrict our discussion to the low-energy molecular 
orbitals which govern the transport phenomenon around 
the Fermi level (Ep). Note that, since we have incor- 
porated a large number of Au atoms to model the ex- 
tended molecule, the relative location of the molecular 
levels with respect to Ep is accurately estimated. On 
examining the T(E) for all the cases in the energy range 
Ep-2 eV to Ep+2 eV, we find that the transmission for 
the non-protonated basepairs is governed by the occupied 
molecular orbitals whereas for the protonated basepairs, 
the transmission is controlled by the unoccupied molec- 
ular orbitals. We note that this is a direct consequence 
of the positive AEA (Adiabatic Electron Affinity) asso- 
ciated with the free basesfTS]. 

Before analyzing the T(E) data for the individual cases 
we point out the general features in the T(E) data plot- 
ted in Fig. 2 . The naturally occurring Watson -Crick 
basepairs A:T and G:C have an inherent asymmetry 
due to which we find that the transmission peaks fol- 
low the changes in the chemical potential of one elec- 
trode. Whereas, for the symmetric basepair A:A, we 
find no such preference and the transmission peaks fol- 
low a symmetric trend in both the forward and reverse 
bias directions. Similar shifts of molecular levels under 
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FIG. 3: (color online) (a) Rectification mechanism for A:Tp. 
(b) Switching action in noncanonical basepair A:A. Shaded 
portion illustrates the integrating region. Arrows indicate the 
bias direction. Solid lines indicate strong transmission and 
dotted line indicate weak transmission. 



the influence of an external bias have been studied be- 
fore for a symmetrically and asymmetrically contacted 
molecule [TrJ]. 

The HOMO band of G:C NP (Fig.2(a)) which is 0.35 
eV below the Ep shows very weak zero bias transmission 
and it does not improve with a change in the bias volt- 
age. However, the HOMO-1 band which is 1.35 eV below 
the Ep, shows strong transmission in the negative bias 
region, but not in the forward bias region. Also we note 
that in the negative bias region, the transmission peak 
moves away from the Ep and thus we do not expect this 
transmission band to enter the integration region and 
have any pronounced effect on the Current- Voltage (I- 
V) characteristics of the system. For G:Cp (Fig. 2(b)), 
we do not find any strong transmission peaks in the en- 
ergy range considered by us. Thus, our calculations show 
that the G:Cnp and G:Cp do not posses any features of 
interest for molecular electronics applications. 

However, the situation is very different for the A:T 
systems. As can be seen from Fig. 2c, the HOMO band 
for A:Tjvp, which is 0.75 eV below the Ep, shows strong 
transmission in the negative bias region, but not in the 
forward bias region. However, due to the asymmetric 
nature of the electron distribution in the basepair, the 
transmission peak moves away from the Ep with increase 
in bias, in the negative bias region. Thus, we do not ex- 
pect this transmission band to enter the integration re- 
gion to generate any nonzero I-V characteristics for small 
reverse bias. However, for the A:Tp (Fig. 2(d)), we ob- 
serve that the LUMO band which is 0.47 eV above the 
Ep shows very strong transmission at zero bias and ap- 
preciable transmission in the negative bias as well as in 
the forward bias region. We also find that in the positive 
bias region, the transmission peak appears right at the 
Ep . 

From the above observations, we describe a rectifica- 
tion mechanism for the A:T basepair (Fig. 3(a)). In the 
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negative bias region, the HOMO level for A:Tnp and 
the LUMO level for A:Tp do not enter into the integrat- 
ing region (illustrated by the shaded region in Fig. 3(a)). 
However, in the positive bias region, low-energy levels, 
namely, HOMO level of A:T NP and the LUMO level of 
A:Tp enter the integration region and thus we would ex- 
pect an increase in the current signal for both A:Tjyp and 
A:Tp for this bias polarity. Interestingly however, due to 
the weak transmission (indicated by the dotted lines in 
Fig. 3a) associated with the HOMO level for A:Tjyp in 
the positive bias region, the resulting current would be 
quite small. Thus, we note that, only in the A:Tp sys- 
tem, we will observe an asymmetric rectifying I-V profile 
with significant current in the positive bias polarity and 
negligible electron flow in the reverse polarity of applied 
bias. 

On analyzing the T(E) data for the A:A systems, we 
find that using the reversible equilibrium between A:Aj\rp 
and A:Ap, the A: A basepair can be used as a molecu- 
lar switch at zero bias. At equilibrium (zero bias), the 
HOMO level of A:Ajvp (Fig.2(e)) is 1.05 eV below the 
E F , whereas the LUMO level of A :A P (Fig. 2(f)) occurs 
right at the Ep. This is quite significant, as even at 
zero bias, if the pH of the system is changed from neu- 
tral to acidic, we would see a switching behavior which 
is reversible. From the T(E) data we also find that this 
switching action can be observed in a very small applied 
bias window as the transmission of the LUMO level of 
AAp decreases rapidly under the influence of an ex- 
ternal bias in both the directions. We depict this sce- 
nario schematically in Fig. 3(b). The mechanism of the 
switching behavior can be explained by examining the 
HOMO and LUMO orbitals of A:Ajvp and A:A P , re- 
spectively. The same are plotted in Fig. 4. We find that 
the HOMO of AA^rp is not delocalized over the entire 
basepair, rather is localized on the individual Adenine 
bases, whereas, the LUMO of A:Ap is delocalized over 
the entire basepair. Note that, such large derealization 
characteristics of A:Ap stems from the 7r-7r interactions 
between the Adenine bases. In a recent experimental 
study by Seferos et al. the strength of tt-tt coupling ver- 
sus conjugated coupling was studied and it was found 
that the tt-tt coupling improved the conduction due to 
strong through space electrostatic interactions. This is 
the same reason why A:Ap shows stronger conducting 
profile in comparison to A:Ajvp[17l [18]. Also, the fact 
that the switching interactions is controlled due to the 
interplay between weak interactions, namely, H-bonding 
and tt-tt interactions, promise technological applications 
with easy tunability. Our findings thus suggest that the 
noncanonical AA systems are prospective candidates for 
switching action under very small pH change to the over- 
all system. 

In conclusion, we find that transport characteristics of 
the A:T and the G:C basepairs vary significantly from 
each other. The H-bonds in A:T are energetically very 
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FIG. 4: (Color online) Orbital plots of the basepairs; HOMO 
of A:Ajvp (left) and LUMO of A:A P (right). 



different from those in G:C; while A:T favors efficient 
transport, G:C does not. We find that under mild acidic 
conditions, the A:T basepairs can be used as efficient 
molecular rectifiers; the action of which can be enhanced 
by tuning the pH of the system. Also we find that the 
non-canonical A:A basepair can be used as an efficient 
molecular switch at zero bias, where the switching be- 
havior is controlled by an interplay between weak inter- 
molecular H-bonding and 7r-7r interactions. 
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